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Abstract 

Introduction: This paper describes a method for automatic analysis of the choroid in 
OCT images of the eye fundus in ophthalmology. The problem of vascular lesions 
occurs e.g. in a large population of patients having diabetes or macular 
degeneration. Their correct diagnosis and quantitative assessment of the treatment 
progress are a critical part of the eye fundus diagnosis. 

Material and method: The study analysed about 1 '000 OCT images acquired using 
SOCT Copernicus (Optopol Tech. SA, Zawiercie, Poland). The proposed algorithm for 
image analysis enabled to analyse the texture of the choroid portion located 
beneath the RPE (Retinal Pigment Epithelium) layer. The analysis was performed 
using the profiled algorithm based on morphological analysis and texture analysis 
and a classifier in the form of decision trees. 

Results: The location of the centres of gravity of individual objects present in the 
image beneath the RPE layer proved to be important in the evaluation of different 
types of images. In addition, the value of the standard deviation and the number of 
objects in a scene were equally important. These features enabled classification of 
three different forms of the choroid that were related to retinal pathology: diabetic 
edema (the classification gave accuracy ACCj =0.73), ischemia of the inner retinal 
layers (ACC 2 = 0.83) and scarring fibro vascular tissue {ACC 3 = 0.69). For the cut 
decision tree the results were as follows: ACC ] = 0.76, ACC 2 = 0.81 , ACQ = 0.68. 

Conclusions: The created decision tree enabled to obtain satisfactory results of the 
classification of three types of choroidal imaging. In addition, it was shown that for 
the assumed characteristics and the developed classifier, the location of B-scan does 
not significantly affect the results. The image analysis method for texture analysis 
presented in the paper confirmed its usefulness in choroid imaging. Currently the 
application is further studied in the Clinical Department of Ophthalmology in the 
District Railway Hospital in Katowice, Medical University of Silesia, Poland. 

Keywords: Eye, Image processing, OCT, Texture analysis, Conditional erosion and 
dilation 



Introduction 

Choroid plays an essential role in many physico-chemical processes. The structure is 
important also for ciliary-retinal vessels (observed in minority of patients) originating 
from the choroid which supply the speckle field and protect against loss of central vi- 
sion, for example in the case of central retinal artery (CRA) occlusion [1]. Visible chor- 
oidal vessels are found to a lesser extent in foveal avascular zone (FAZ). In the case of 
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fluorescein angiography, it shows no presence of the fluorescence in that area (due to 
high amount of pigment). Conversely, any leakage of the pigment or FAZ staining indi- 
cates macular disease. Fluorescein diagnostic method is, in this case, profiled to carry 
out this type of diagnosis [2]. In practice, therefore, diagnosis of the eye fundus and 
choroidal layer using optical coherence tomography (OCT) also brings correct results 
[3,4]. The analysis of the vascular layer located beneath the RPE layer (retinal pigment 
ephitelium) in the OCT image of the eye is presented in a small number of publications 
[4-14]. They are mainly related to qualitative analysis of the choroidal layer without 
quantitative treatment of the characteristic distributions which occur there. Therefore, 
in this paper quantitative analysis of the choroidal layer using the new developed algo- 
rithm for image analysis and processing was proposed. Most of currently used OCT de- 
vices are not intended for choroidal imaging. The authors tried to obtain data resulting 
from the choroid reflectivity in a wide variety of patients. The algorithm was profiled to 
the analysis of the following types of images: 

- neovascular AMD or exudations secondary to diabetic or thrombotic edema 
(specific layouts of shadows in the choroid caused by retinal changes) [5,6], 

- diffuse macular edema without blood and exudations or ischemia of the inner 
retinal layers (global reduction of brightness in an OCT image) in such patients 
there is need to differentiate between the choroidal atrophy due to degeneration or 
high myopia [7], 

- scarring fibrovascular tissue a uniform image proves its presence in most patients [7-9]. 
These types of images are presented in Table 1 they are subject to further analysis. 

Material 

The study analysed about l'OOO OCT images acquired using SOCT Copernicus 
(Optopol Tech. SA, Zawiercie, Poland). The patients ranged in age from 12 to 78 years 
and had different types of choroidal structure. It was a group of patients routinely exam- 
ined, analyzed retrospectively and anonymously. The routine tests were carried out in 
accordance with the Declaration of Helsinki. The images were acquired in DICOM or 
RAW format with a resolution of 256x1024 pixels at 8 bits per pixel. Image analysis 
was carried out in Matlab with Image Acquisition Toolbox and Signal Processing tools 
(version 4.0 and 7.1 respectively), whereas code optimization was carried out in the C 
language. The proposed algorithm for image analysis enabled to analyse the texture of 
the choroid portion located beneath the RPE (Retinal Pigment Epithelium) layer. The 
analysis of the choroid was performed using the new profiled algorithm based on tex- 
ture analysis and mathematical morphology that is described below. The division into 



Table 1 Types of images and their features visible in OCT images 



Symbol 


Imaging type 


Features visible in the image 


z, 


neovascular AMD or exudations secondary to diabetic 


characteristic layouts of shadows in the 




or thrombotic edema 


choroid caused by retinal changes 


z 2 


diffuse macular edema without blood and exudations 


global reduction in brightness in the OCT 




or ischemia of the inner retinal layers 


image 


z 3 


scarring fibrovascular tissue 


uniform image is evidence of its presence 
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types of images (3 groups) was carried out by an ophthalmologist for 1 000 images 
representing learning, validation and test groups (proportion: 60%, 20%, 20% respect- 
ively). At this stage, patients with other type of images were also eliminated from fur- 
ther analysis - only images with visible alterations in the choroid were further analyzed. 

Method 

Preprocessing 

The input image L GRAY with a resolution MxN = 256x1024 pixels (M - number of 
rows, N - number of columns of the image) was subjected to median filtering using a 
mask h sized M h xN h = 3x3 pixels, thus obtaining the image L MED . The size of the filter 
mask h was chosen on the basis of medical evidence on the extent of artefacts found in 
this layer of the eye fundus OCT image and resolution of the image L GRAY . The images 
of successive stages of image pre-processing are shown in Figures 1 and 2. The image 
L MED is further analysed to detect the RPE layer Jrp E (h) [3,15-17]. For this purpose, 
every «-th column of the image L MED is examined. The position of maximum bright- 
ness for each column is determined, i.e.: 

yRPE{ n ) = ar S max L MED (m,n) (1) 

me(l,M) 

where: 

m,n - coordinates of rows and columns of the matrix ms(l,M) and ne(l,N). 

The Equation (1) can be directly applied only if for all the analysed rows and succes- 
sive columns, there is only one maximum value of brightness. In practice, it occurs in 
about 80% of the analysed cases at the resolution of 8 bits per pixel. The Equation (1) is 
also very sensitive to noise, especially in the case of single bright pixels - salt and pep- 
per noise. This type of noise is sometimes not fully filtered during median filtering. For 
this reason, the Equation (1) was expanded to the following form (y R p E " («)): 

(m if L MED (m,n) < max L MED (m,n)*p r 
L MA x(m,n)=< me(iM) (2) 

[ 0 others 

{M 
med (LMAx(m,n)) if VZmax(w,«) _ 
me(l,M)m*0 ^ > U W 

M other 

where: 
med - median 

p r - coefficient determined in the range from 0 tol. 

The value of the coefficient p r [18] is determined once for the Equation (2) and is 0.9 
(it was determined arbitrarily). Depending on its selection (p r ), the number of pixels 
that influence the calculation of the median is changed. An increase in the value of p r 
increases the number of pixels from which the median is calculated. Equations (2) and 
(3) by selecting p r enable correct calculations of the RPE layer individual points only if 
there is one cluster. In other cases, when there are two or more clusters, the calcula- 
tions are more difficult. These are cases where the RPE layer is not the brightest layer for 
the analysed column. Such situations are very rare. The specificity of the Equations (2) 
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RPE 



CHO 



Figure 1 The method for obtaining tomographic images of the fundus. For a sample 2D tomographic 
image, the RPE layer (retina pigment epithelium) and the choroid layer CHO are highlighted. Image analysis 
applies to the proposed algorithm which analyses the choroid layer using new methods of texture analysis 
and mathematical morphology. In each case, a flat two-dimensional input image is analysed, whose reso- 
lution (and that of the OCT apparatus) does not affect the obtained results. 



and (3) enables to receive values equal to M (last row) for cases where none of the pixels 
is larger than max m E (ij^L MED (m, n) * p r for the analysed column. 

Based on the course oi y RPE " {n) (hereinafter abbreviated to y RPE ), the image L RPE 
was created in the following way: 



z(m,n) 



(L M ED(m,n) if m>y RPE (n) 
[ -1 others 



(4) 



The values of added to the matrix L RPE enable to distinguish the area above the 
RPE layer from the pixels beneath the RPE with a value of "0". The final element of 
image pre-processing is affine transformation of the image L RPE to the image L CHO 
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Figure 2 The course of image pre-processing. The subsequent steps of L GRAY image analysis: filtration 
with a median filter - L MED , determination of the RPE layer (retinal pigment epithelium) - ynpE(n), coordinate 
system conversion - L C ho- These steps are part of the image pre-processing which is necessary for further 
analysis of images. 



containing only the interesting area of the choroid Figure 2. The image L C ho includes 
all the pixels of the image L RPE from the layer y R pE(n) to the last row (M). The 
remaining space in the matrix L RPE is filled with the values Thus prepared image 
L C ho with a resolution M c xN c is subjected to appropriate analysis and processing as 
described in the next sub-section. 

Image processing 

The input image L CHO shown in Figure 2 provides the basis for further processing. The 
specific properties of optical scanners as well as the object (eye) specificity contribute 
to the fact that brightness uniformity correction is necessary for further processing. For 
this purpose, the image L MEAN , which results from filtration with an averaging filter 
with a mask h 2 sized M h2 xN h2 = 30x30 pixels, was subtracted from the image L CHQ . As 
a result, the image L CHOM was obtained, i.e.: 

L C HOM(m, n) = L C Ho(m, n) 

1 srsrfr ( , m„ 2 N h2 \ \ 

-^LL {Lena [m + m 2 - — ,n + n 2 -—j.h 2 (m 2 ,n 2 )j 

(5) 

for me(M h2 /2, M c -M h2 /2) and ne(N h2 /2, N c -N h2 /2). 

The size of the mask h 2 is closely related to the size of objects subjected to detection. 
In the case of the diseases listed in Table 1, the maximum size of objects is 20x20 
pixels. Therefore, it was assumed that the size of the mask h 2 must be at least twice as 
large. This is necessary to carry out the removal of uneven brightness. Further process- 
ing steps are related to the detection of objects with the use of morphological opera- 
tions. For these operations, it is necessary to define the binary image L BW . This image 
contains information on the location of pixels of the image L CHOM derived from the 
image L RPE i.e.: 
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i \ f 1 if L RPE (m,n) = -1 

={ 0o i rs (6) 

The proper steps of image processing are related to sequential morphological analysis. 
Morphological opening operations were performed on the image L CHOM . The size of the 
structural element S£, was varied every 2 pixels in the range from 3x3 to 11x11 pixels 
(where the subscript i indicates the size, i.e.: SE 3 is a structural element sized 3x3 etc.). An 
opening operation was carried out for every z'-th size of the symmetric structural element 
SE, i.e.: 

L 0i = min max{L C HOM) I -L BW (7) 

SEj \ sEj J 



These images (L Ql ) after normalization to the range 0 to 1 may be subjected to binari- 
zation at a constant threshold p r . 



1 if I L 0i - min(I 0i ) / max I L oi - min (L oi ] 

\ Loi J Loi \ La 

0 others 



< Pr (8) 



The value of the threshold p r is fixed at 0.5 due to the previously performed opera- 
tions of removing uneven brightness and standardization. The images L Bi require ad- 
justments related to the removal of small artefacts and holes in objects. This process 
was carried out using relationships of conditional erosion and dilation of the binary 
image L Bi . In the case of a symmetrical structural element SE it relationships of condi- 
tional erosion and dilation [19] are simplified to the following form: 

L E(c) (m,n) = 

(L Bi (m,n) for (\-p we )-p mn (m.n)<s re {m,n) (g) 

- i min (L Bi (m + m SEh n + n SEi )) for (l-p we )-p mn (m,n) > s re {m,n) 

\m SEi ,n SB ^SEi 



L D(c) {m,n) = 

(L Bi {m,n) for (p wd + l)-p mn (m,n)>s rd (m,n) ( 10 ) 

- i max (L Bi (m-m SE i,n-ns E i)) for (p wd + \)-p mn {m,n) < s rd (m,n) 

y msEi,nsEi^SEi 

where: 

LE(c)( m > n ) - the resulting binary image after subjecting the image L Bi to conditional 
erosion, 

LD(c)( m >n) - the resulting binary image after subjecting the image L B i to conditional 
dilation, 

p we - constant erosion effectiveness, 

p wd - constant dilation effectiveness, 

p mn (m,n) - threshold dependent on the coordinates m, n, 

s re - the mean value of the analysed area for erosion, 

s rd - the mean value of the analysed area for dilation. 

The mean values s re , s rd for erosion and dilation respectively were calculated from the 
following equations: 
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MsEi NsEi J ( 

s re {m,n) = 2^ 

m s£i =l K S£ ,=1 



m SE i,n + n S Ei) 



MsEi-N SE i 



(11) 



v v^t L 0 i(m-msEi,n-nsEi) 
s rd (m,n)= ^ ^ 



(12) 



The constants and determining the effectiveness of erosion and dilation re- 
spectively take the following values: 

/V e e{-1.0, -0.9, -0.1, 0.0, 0.1, 0.9, 1.0}, (13) 
PvM-1.0, -0.9, -0.1, 0.0, 0.1, 0.9, 1.0} 

The choice of this range results from the condition of the left side of inequality, i.e.: 




/ ii 



The values of p mn are in the range from 0 to 1, whereas the values of (l-p we ) should 
be non-negative in the range from 0 to 2. The values of {l-p we ) and (1 + p wd ) for p we = 
p w d = 0 are equal to 1, which means high intensity of conditional operations. For the 
other values of the thresholds p we and p wd , for example for p we = p wd = 1, there is a 
complete lack of effectiveness of erosion operations and significant effectiveness of dila- 
tion (the impact of selection of p we and p wc i values is presented later in this section). 
However, very often p mn {m,n) = const, irrespective of the location (p mn *f[m,n)). Adopt- 
ing p mn {m,ri) = const is due to the nature of conditional operations, where in a general 
case a condition may not only be dependent on the mean values of s re and s,^, but also 
on other values of the pixel saturation degree. These special properties of conditional 
dilation and erosion enable to obtain effective correction of the quality of the input im- 
ages L Bi . Sequential execution of conditional dilation and erosion (in this case, three 
times) allows to obtain corrected images L Ki (Figure 3). The shape of the structural 
element SEi adopted in all of these relationships was as a circle of a pre-specified size 
because of the shape of the recognized objects. For each image L Ki , characteristics were 
determined for each object. These features include: 

w(l) to w(5) - number of objects in the image L Ki for f£(3, 5, 7, 9, 11), 

m/(6) to vw(10) - the average position of the centre of gravity in the x-axis for all the 

objects in the image L Ki for z'£(3, 5, 7, 3, 11), 

vv(ll) to w(15) - the average position of the centre of gravity in the j-axis for all the 
objects in the image L Ki for z'£(3, 5, 7, 9, 11), 

w(16) to w(20) - standard deviation of the mean brightness of pixels of all the objects 
in the image L Ki for /£(3, 5, 7, 9, 11), 

These features were selected taking into account medical conditions. They are related 
to the location of vascular lesions from w(6) to w(15), uniformity of brightness distribu- 
tion within these changes from w(l6) to w(20) and the area of changes for the appro- 
priate number of objects from w(l) to w(5). These characteristics form the basis for 
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Figure 3 Image analysis of the sequence of images L Ki . Subsequent results are shown for /e{3, 5, 7, 9, 
1 1). For each image L KI and thus for each / the values of the features from w(1) to w(20) are calculated. For 
example, one of the objects whose coordinates of the centre of gravity are calculated is shown on the top 
of the zoom - in this case, they are (1 76, 281). The features w(6) to w(16) are the mean value of gravity 
centre coordinates of all objects in the image L Ki . 



building a classifier using decision trees (Figure 4 shows the block diagram of the 
algorithm). 

Results 

The values of the 20 features obtained (four different types) are further used to build a 
decision tree. An ophthalmologist divided more than 1 000 images into three groups 
(Z lt Z 2 or Z 3 - Table 1). The division into groups Z h Z 2 or Z 3 was performed manually 
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o 



Decision tree cut 



Classification Zi, Z2, Z3 



Figure 4 The block diagram of the algorithm. The block diagram is divided into three main parts, 

namely image pre-processing using the known techniques of image analysis and processing, the principa 

analysis of images proposed by the authors and the classifier based on decision trees. The individual 

processing blocks are described in detail in the paper, 
s * 
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by an ophthalmologist. As a result, the following frequency of occurrence of choroidal 
imaging was obtained: the group Zj included 23% of all cases, the group Z 2 included 
44% of all cases, and the group Z 3 included 33% of all cases. This corresponds to 233, 
436 and 331 images respectively. Whereas the division into learning, validation and test 
groups had the following proportions: 60%, 20% and 20% respectively. These groups 
were formed after rejecting the images with a mixed character of the changes observed - 
the type of the disease was not a criterion for exclusion here. A variety of overlapping 
diseases can be found in the excluded images. The cases of spatially invisible layer of 
choroid, images resulting from the errors in the acquisition or lacking a visible layer of 
RPE (for various reasons) were also excluded. Due to the fact that these are retrospect- 
ive studies, the excluded images did not often cover the full range of the choroid, were 
deliberately obscured or distorted at the acquisition stage. 

The nodes of the decision tree are different features from w(l) to w(20), the branches 
are the values corresponding to these attributes, and the leaves make individual deci- 
sions - type identification {Z lt Z 2 or Z 3 - Table 1). In all cases, a non-parametrical al- 
gorithm creating CART (Classification and Regression Trees) binary trees was used as 
the method for decision tree induction. An increase in the nodes purity was used as 
the criterion assessing the quality of CART divisions. The Gini index was used as the 
measure of nodes impurity. The tree creation was not limited by a minimum number 
of vectors in a node. As the considerations apply to the construction of a classifier 
based on the knowledge base, w(l) to w(20) features, a preliminary prepared tree 
Figure 5 was built based on the full information, using the training group. The results of 
the classification for the complete decision tree for the test group as follows: 

- for Zj - SPCj = 0.74, TPRj = 0.32, ACCj = 0.64, 

- for Z 2 - SPC 2 = 0.56, TPR 2 = 0.86, ACC 2 = 0.70, 

- for Z 3 - SPC 3 = 0.88, TPR 3 = 0.097, ACC 3 = 0.63. 

(where: SPC = TN/(FP + TM) - specificity, TPR = TP I {TP + FN) - sensitivity, TN - 
true negative, TP - true positive, FN - false negative, FN - false positive, subscripts 
"1", "2" or "3" indicate the type of images). A closer analysis of the resulting decision 
tree enables, at this stage, rough assessment of the importance of individual features. 
For this form of the complete decision tree, there occur only features w(l), w(6), w(15) 
and w(20). Additional information is provided by the ROC graph (Receiver Operating 
Characteristic) designated for Z lt Z 2 , Z 3 to assess the impact of individual features from 
w(l) to w(20) separately. This graph (ROC) is shown in Figure 6. Based on the ROC 
graph and the complete decision tree, it can be concluded that there is a negligible im- 
pact of the features w(ll) to w(20) on the obtained results. Additionally, the sensitivity 
of features to changes in the decision-making threshold value is the smallest for the fea- 
ture w(ll), and the highest for w(8) (Figure 6). For this reason, decision trees were created 
for each configuration of features, thereby forming 2 20 full decision trees (unpruned trees) 
and the same number of cut decision trees (over two million decision trees in total - 2' 
097 '152). Figure 7 shows the graph of cost relationships (misclassification error) as a func- 
tion of stratified cross-validation and resubstitution. The graph was also shown by com- 
puting a cutoff value that is equal to the minimum cost plus one standard error. The best 
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Figure 5 The complete classification tree. The decision tree was created for the classification of imaging 
types Z h Z 2 , Z 3 for the features to w(20). The characteristic elements such as the presence of the feature w 
(20) in the first node of the tree and successively features w{6) and w(15) for the other nodes are visible 
here. The results obtained on this basis are, for example, SPC 2 = 0.56, TPR 2 = 0.86, ACQ = 0.70 for Z 2 . 



level computed by the classregtree test method is the smallest tree under this cutoff (best 
level = 0 corresponds to the unpruned tree). The best results obtained for the complete 
decision tree are shown in Table 2. The best results of ACC 2 = 0.83 were obtained for the 
features w(l), w(2), m/(3) and w{4). 

In the classification of groups Z lt Z 2 , Z 3 the features from w(9) to w(12) as well as from 
w(l8) to w(20) did not occur. This means that their influence can be neglected for the best 
classifications (top 3 results). These features (from w{9) to tv(12) as well as from w(18) to 
w(20)) define the average location of the centre of gravity for i = 9 and 11 in the x-axis, for 
i = 3, 5 in the j-axis and the standard deviation of the average brightness of pixels of all 
objects for i = 9 and 11. The same results of ACC 2 were obtained for various configura- 
tions of the features w(l), w(2), w(3) and w(4?) or w(5), w{7), w(8), w(16) and w(17), which 
is also quite interesting (Table 2). The best results for the cut decision tree are shown in 
Table 3. ACC 2 = 0.82, obtained for the features w(l), w(4) w(5), w(7) w(8), w(12) and w 
(15), is the best result for the cut decision tree. The results for the groups Z lt Z 2 , Z 3 con- 
firm that the features from w(17) to w(20) have the smallest influence. The summary fre- 
quency chart of the occurrence of individual features from w(l) to w(20) for the cut and 
complete decision trees for the first l'OOO best results is interesting (Figure 8). It shows 
the greatest frequency of occurrence of the feature w(14) for the complete decision tree 
and the feature w(4) for the cut decision tree. From a practical standpoint, however, min- 
imizing the number of features occurring in the classification seems to be the most vital. 
Therefore, when analysing only the three top results for the cut decision tree, one feature 
for Zj, six features Z 2 and one feature for Z 3 are obtained. This means that only the group 
Z 2 requires the largest number of features in the classification. This information is essen- 
tial for optimizing the computational complexity of the algorithm. The analysis time for a 
single image does not exceed one second for the Pentium 4 CPU 3.0 GHz, 8GB RAM. 
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1-SPC 

Figure 6 ROC Chart (Receiver Operating Characteristic). The graphs show the impact of threshold values 
for individual features w(1) to w(20) for the different imaging types Z h Z 2 , Z 3 . A negligible impact of the features 
w(1 1) and w(20) on the obtained results can be inferred here. In addition, sensitivity of features to changes in 
decision-making threshold values is the smallest for the feature w(] 1 ) and the highest for w(8). 



Comparison with other authors' results 

The analysis of the choroidal layer is dealt with in a number of publications [20-23]. In 
most cases, however, it is qualitative analysis. It includes the analysis of the choroidal 
layer thickness [20] or analysis with the use of Heidelberg Eye Explorer software [21]. 
In the first case, central serous chorioretinopathy case-series is presented. The follow- 
ing results were obtained - mean subfoveal choroidal thickness was significantly (sig- 
nificance level p s i = 0.04) larger in the affected eyes (455 ± 73 um) than in the 
contralateral unaffected eyes (387 ± 94 um), in which it was significantly (p s! = 0.005) 
larger than in the control group (289 ± 71 um). The control group and the group of pa- 
tients in this case [20] comprised 15 + 15 patients. The largest vessel diameter was 
significantly (p s i < 0.001, correlation coefficient: 0.73) correlated with the thickness of 
the total choroid. A similar small number of patients was analysed in [21]. The de- 
scribed software enables to measure some features of the choroid. However, this 
method is fully manual. Another interesting analysis is shown by Mrejen S. in [22]. The 
results obtained there concern the use of various imaging techniques in the diagnosis 
of the choroid. An automatic method for the analysis of the choroid is shown in the 
work of Park S.Y. [23]. The analysis presented there concerned only determination of 
the location of the selected layers in the tomographic image of the fundus. In addition, 
the correlation between the early treatment diabetic retinopathy study and EDI (en- 
hanced depth imaging) was demonstrated. The automatically measured retinal thick- 
ness and volume of 9 early treatment diabetic retinopathy study subfields with 
conventional and EDI raster scan showed an intraclass correlation coefficient of 0.861 
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Figure 7 The graph of the cost relationships (misclassification error) as a function of stratified 
cross-validation and resubstitution. The blue shows the cross-validation error. The best choice, the place 
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error (min + 1 std. err). 



Table 2 The three top results obtained for different configurations of the features w(1) 
to w(20)* for the complete decision tree 



w/N° 




z, 






z 2 






z 3 




1 


2 


3 




2 


3 


1 


2 


3 


w(1) 


0 


0 


0 






0 


1 


1 


1 


M2) 


0 


0 


0 






0 


0 


0 


1 


M3) 


0 


0 


0 






0 


1 


0 


0 


M4) 


I 


1 


1 






0 


0 


0 


0 


MS) 


0 


0 


0 


0 




1 


0 


0 


0 


w(6) 


1 


1 


1 


0 


0 


0 


1 


0 


0 


w{7) 


0 


0 


0 


0 


0 


1 


0 


1 


1 


w{8) 


0 


1 


0 


0 


0 


1 


0 


0 


0 


M9) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(10) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(11) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(12) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(13) 


0 


0 


1 


0 


0 


0 


0 


0 


0 


w(14) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(15) 


1 


1 


1 


0 


0 


0 


0 


0 


0 


w(16) 


0 


0 


0 


0 


0 


1 


0 


0 


0 


w(17) 


0 


0 


0 


0 


0 


1 


0 


0 


0 


w(18) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(19) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(20) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


ACQ 


0.73 


0.73 


0.73 


0.83 


0.83 


0.83 


0.69 


0.69 


0.69 



*The occurrence of a given feature w is indicated by "1" and its lack by "0". 
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Table 3 The three top results obtained for different configurations of the features w(1) 
to w(20) for the cut decision tree 



w/N° 










z 2 






Z 3 




1 


2 


3 


1 


2 


3 


1 


2 


3 


w(D 


0 


1 


0 


1 


1 


1 


1 


0 


0 


w(2) 


0 


0 


1 


0 


0 


0 


1 


0 


0 


w(3) 


0 


0 


0 


0 


1 


1 


0 


1 


0 


w(4) 


I 


1 


1 


1 


0 


0 


0 


0 


1 


w{S) 


0 


0 


0 


1 


0 


1 


0 


0 


0 


w(6) 


0 


0 


0 


0 


1 


1 


0 


0 


0 


w(7) 


0 


0 


0 


1 


0 


0 


0 


0 


0 


w{8) 


0 


0 


0 


1 


0 


0 


0 


0 


0 


w(9) 


0 


0 


0 


0 


1 


1 


0 


0 


0 


w(10) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(11) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(12) 


0 


0 


0 


1 


1 


1 


0 


0 


0 


w(13) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(14) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(15) 


0 


0 


0 


1 


0 


0 


0 


0 


0 


w(16) 


0 


0 


0 


1 


1 


1 


0 


0 


0 


w(17) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(18) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(19) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


w(20) 


0 


0 


0 


0 


0 


0 


0 


0 


0 


>4CC, 


0.76 


0.76 


0.76 


0.82 


0.81 


0.81 


0.68 


0.68 


0.68 



to 0.995 and 0.873 to 0.99 respectively. The number of analyzed cases was 35 patients 
with chorioretinal diseases and 20 healthy subjects. In this paper, the number of cases 
is much higher and the image analysis and processing are fully automatic. It seems that 
useful results can be obtained even with non-EDI SOCT. However, it should be further 
tested with other devices. 

Summary 

This paper shows a new fully automated method for the analysis of choroidal images. A div- 
ision into three types of choroidal images was proposed. The created decision tree enabled 
to obtain satisfactory classification results. For example, for the classification of images Z 2 , 
the accuracy ACC 2 was 0.81 (pruned tree). Additionally, the features which have the great- 
est impact on the classification efficiency were characterized (creating more than 2 million 
decision trees). These are the characteristics from w(l) to w(9), w(12), w(15) and w 
(16) which are responsible for the number of objects in the scene and the average pos- 
ition of the centre of gravity in the x-axis for i = 3, 5, 7, 9 and the average position of the 
centre of gravity in the y-axis for i = 5, 11 and STD for i = 3. In addition, it was shown that 
for the adopted classifier of the cut decision tree, for the top 1 000 results the feature w 
(4) was most common. It was also noted that the different locations of scanning 
differentiate the images of the choroid but the classification result remains the same. 
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Figure 8 The graph of occurrence frequency of individual features w(1 ) to w(20) for the complete 
decision tree and the cut one for the first top 1'000 results. The graph shows that the highest occurrence 
concerns the feature w(14) for the complete decision tree and the feature w(4) for the cut decision tree. The 
occurrence frequency of the features corresponds with the results shown in Tables 2 and 3. 

Foveal location, however, has the greatest ability to differentiate changes. Individual 
B-scans located on the macula periphery are closer to each other. The proposed al- 
gorithm for image analysis and the classifier were implemented in the C language. Cur- 
rently, the application is further studied in the Clinical Department of Ophthalmology in 
District Railway Hospital, Medical University of Silesia in Katowice, Poland. 

The presented methodology of the procedure does not cover a large range of opportun- 
ities offered by modern methods of image analysis and processing. Currently, there is on- 
going work on the optimization of the algorithm in terms of computational complexity. 
Moreover, its aim is to replace some of the time-consuming functions with some other 
more efficient ones. In particular, the possibility of applying different methods of texture 
analysis [23-31], Boolean function [32] or other methods of image analysis [33,34] or the 
impact of other factors [35-41] is considered. Other types of classifiers ensuring the ana- 
lysis of other characteristics acquired from images are particularly taken into account. 

Consent 

Written informed consent was obtained from the patient for the publication of this 
report and any accompanying images. 
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